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SUMMARY 


The  propagation  of  continuous  plane  progressive  sound  waves 
whose  pressure  variation  is  of  the  order  of  one  tenth  of  the  average 
pressure  is  discussed.  It  Is  shown  that  regardless  of  the  initial  wave 
form  shocks  will  develop  at  the  leading  front  of  each  wave  after  sev- 
eral wave  lengths  of  propagation.  Assuming  the  resulting  stable  wave 
form  to  be  saw-tooth  in  character  the  attenuation  of  these  repeated 


shocks  is  derived  from  shock  w&Vh  theory.  Writing  Rj/p^  ° 1 , 

where  p — p is  the  pressure  discontinuity  at  the  shock,  it  is  shown 

u 1 


that 


I _ •/**  I X ~Xe  . 

" To  - Tx  • “T 


where  ^ is  the  value  of  d at  the  distance,  \ } equal  to,  , / 

is  the  ratio  of  specific  heats  and  \ is  the  wave  length  of  the  sound. 
This  result  is  discussed  and  compared  with  previously  published  studies 
of  tho  attenuation  of  single  N waves  and  found  to  be  compatible.  Also 
it  is  shown  that  Fay's  solution  of  the  hydrodynamic  equations  including 
the  effects  of  viscosity,  which  shows  the  stable  wave  form  to  be  a saw- 
tooth, may  be  extended  to  yield  the  attenuation  rate  derived  here. 
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The  Attenuation  of  Very  High  Amplitude  Sound  Waves 
L Introduction 
The  problem  of  concern  here  1b  the  attenuation  experienced 
by  a very  high  amplitude  progressive,  plane  sound  wave  — on©  whose 
peak  to  trough  pressure  variation  la  of  the  order  of  one  tenth  of  the 
mean  pressure.  Consider  then  a wave  whose  pressure  variation  Ini- 
tially can  be  described  by  a sine  function.  As  la  well  known*,  the 
wave  will  be  distorted  as  it  propagates  since  the  peaks  travel  at  a 
higher  velocity  than  the  troughs  and  consequently  there  will  be  a ten- 
1 dency  for  the  gradient  to  Increase  on  the  leading  front  of  each  wave 
and  decrease  on  the  trailing  front.  Fay2  has  shown  that  a limiting 
form  of  the  wave  when  one  Includes  viscosity  as  an  attenuative  process 
Is  saw-tooth  with  the  leading  front  approaching  a vertical  slope.  It 
would  be  possible,  by  means  of  an  Integration  uf  the  hydrodynamic 

3 

equations  , to  trace  out  In  detail  the  process  of  this  deformation  up 
to  the  development  of  the  shock  front,  but,  for  the  purposes  of  our 
discussion,  all  that  is  needed  Is  a treatment  which  will  yield  the  ap- 
proximate distance  of  propagation  necessary  for  this  deformation  to 
become  complete.  To  obtain  this  it  is  necessary  to  recognize  that 
the  instantaneous  velocity  of  propagation  of  a given  part  of  the  wave  Is 
equal  to  the  sum  of  the  sound  velocity  at  the  given  phase  of  the  wave, 


*See,  e.  g.,  pp.  30-41  of  Ref.  1. 
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c,  and  the  local  particle  velocity,  u.  Thus,  In  Figure*-  (1), 


«tire.choA  of 


p^opa^ct 


\ i0»\ 


Ft  i. 


which  represents  a sine  wave,  the  point  marked  A will  travel  faster 
than  3 because  the  temperature  at  A is  elevated  and  consequently 
c^>  c^  , and  because  the  particle  velocity,  u,  in  tha  direction  of  prop 
agation  is  greater  than  that  at  B.  For  an  adiabatic  wave, 

£ + Li  s.  £g  f — — ' t 

where  i is  the  ratio  of  specific  heats. 

A sufficient  approximation  (for  the  present  purposes)  will  be 
to  assume  that  the  relation  between  jp1  and  u is  that  which  obtains  for 
infinitesimal  sound  waves,  and  that  u at  the  peak  remains  constant  as 
it  distorts.  V/ith  this  approximation  it  is  possible  to  calculate  the 
distance  of  propagation  necessary  before  point  A overtakes  point  B, 

If  r is  the  sound  velocity  at  point  E,  the  density  of  gas  at  point 
B,  then  this  distance,  X,  will  be  given  by 
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It  Ls  thus  clear  that  for  excess  pressure  amplitudes,  p', equal 
to  0, 05  of  tho  ambient  atmospheric  pressure,  it  is  reasonable  to  expect 
shock  wave  character,  at  the  very  most,  several  wave-lengths  from 
the  source, , and,  moreover,  this  result  is  correct  as  to  order  of  mag- 
nitude regardless  of  the  form  of  the  initially  produced  wave.  Conse- 
quently, of  prime  importance  in  considering  the  attenuation  of  very 
high  amplitude  sound  v/aves  is  the  attenuation  of  repeated  shock  waves, 
and  it  is  to  this  problem  which  attention  will  bo  restricted. 

There  are  two  ways  of  approaching  the  problem.  The  more 

p 

fundamental  approach  is  typified  by  the  work  of  Fay  , and  involves  a 
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solution  of  the  hydrodynamical  equations  including  the  attenuative  pro- 
cess of  viscosity  which,  as  pointed  out,  limits  the  maximum  pressure 
gradient.  (Without  such  attenuation  one  la  led  to  multiple  valued 
solutions  after  the  wave  crest  overtakes  the  trough. ) In  this  fashion 
Fay  was  able  to  show  that  a "stable"  wave  form  la  approached  which 
Is  saw-tooth.  A second  procedure  (one  developed  In  detail  In  this 
paper  ) which  has  proved  to  be  lucrative,  can  be  used  and  Is  charac- 
terized by  the  application  of  the  Rankine-Hugoniot  shock  relations. 

One  starts  out  by  assuming  the  wave  has  a saw-tooth  character  and 
the  leading  edge  Is  characterized  by  a discontinuity  In  pressure,  den- 
sity, particle  velocity  and  temperature.  Application  of  the  shock  re- 
lations leads  to  the  result  that  there  Is  an  entropy  increase  across 
the  discontinuity,  or  shock,  and  this  entropy  increase  can  be  inter- 
preted as  a space  rate  of  decrease  of  amplitude.  It  Is  to  be  noted 
that  this  procedure  leads  to  a calculation  of  attenuation  rates  without 
specifically  mentioning  any  attenuative  processes.  (The  substitute  for 
an  explicit  formulation  of  these  processes  is  contained  in  the  assump- 
tion of  the  presence  and  stability  of  the  discontinuity.)  It  will  be 
shown  In  a later  section  of  this  paper  that  the  attenuation  obtained  by 
application  of  the  shock  relations  may  also  be  obtained  from  a result 

O 

given  by  Fay4-. 

The  development  given  here  will  be  a first  order  theory  in 
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the  sense  that  all  but  the  lowest  orders  of  significant  terms  will  be 
dropped.  Also  the  waves,  It  will  be  assigned,  are  exclusively  plana 
progressive  waves. 

n.  Theory 

A.  Propagation  of  a Discontinuity. 

The  results  which  will  be  obtained  in  this  section  A can  be 
found  in  many  texts  (see  e.  g.,  Refs.  (3)  and  (5))  but,  for  the  sake 
of  completeness,  a brief  outline  of  the  procedure  will  be  given. 

Consider  a plane  discontinuity  moving  to  the  right  with  a 
velocity  U.  The  subscript  2 refers  to  the  medium  to  the  left  of  the 
discontinuity  and  1 to  the  medium  to  the  right  of  the  discontinuity. 

Pn  , p «=  pressure 

<J  1 

e density 

Sg,  = specific  entropy 

u , u^  ■ particle  velocity 

Tp,  Tj  = temperature 

% , e1  = specific  Internal  energy' 


It  will  be  assumed  that  the  gas  through  which  the  discon- 


tinuity moves  Is  a polytropic  gas,  i.  e.,  p/>  = RT  and  e = C T,  where 

( v 

C is  the  specific  heat  at  constant  volume, 
v 


Suppose  a velocity  of  flow,  — U,  is  superposed  on  the  whole 
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gas.  Then  the  discontinuity  will  remain  fixed;  gas  will  enter  the  dis- 
continuity from  the  right  with  a velocity  xr,  and  leave  with  a 

velocity  Vj,  s U - U.I.  . 


The  three  shock  relations  assert  that  (1)  mass,  (2)  momen- 
tum, and  (3)  energy  are  conserved  for  the  motion.  They  lead  respec- 
tively to  the  three  relations, 


f.v,  » r,  ^ (i) 

p,  + f,  4 1 = /-'»■  t ^ (2) 

+J-v,'-+e.,  p- + 1 (3) 

/?  i n 

\rt  <£,  j €j,  can  be  eliminated  from  these  equations  in 
the  following  way:  Using  the  polytropic  character  of  the  gas  and  the 

definition  of  V m / f R/dy  > fiq.  (3)  can  be  written  in  the 


From  (1)  and  ("), 


/-/  , 

/;  _ /v  /•  A 

^ V ^ A 
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Now  for  a polytropic  gas, 


-c- 


iifjL 

rJ  ti  ' (?) 


Writs 


where  J is  of  the  order  of  0.1. 

Substituting  (6)  and  (8)  into  (7)  and  expanding  it  in  a power 

series,  it  is  found  that  the  lowest  order  non-zero  term  in  is 

, + 

of  the  third  order  in  £ and  is  given  by 

f c p 'Lll~  a ^ 

S - S,  i K.  ,1  yi  0 ■ (9) 


B.  Attenuation  of  a Repeated  Shock  Wave. 

Consider  now  a vibrating  ristoa  in  a semi-infinite  tube  (see 
Figure  (3)).  Assume  it  to  be  vibrating  for  a sufficiently  long  time  so 


P'3-3. 

th-’t  the  variation  of  hydrodyr.amical  and  thermodynamical  quantities 
at  any;  given  point  in  the  tube  is  truly  periodic  and  that  for  all  points 
the  average  mass  flow  is  nero  and  the  average  temperature  is  uniform. 
Then,  if  the  amplitude  of  the  piston  is  such  as  to  generate*  waves  whose 


See,  e.  g. , p.  41  of  Ref.  5. 
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pressure  variation  is  o i the  order  of  one-tenth  of  the  averago  pressure, 
after  traveling  a few  wave  lengths  from  the  piston  the  fronts  will  have 
developed  shocks,  Assume  that  the  pressure  wave  form  at  this  dis- 
tance is  saw-tooth  and  can  be  pictured  at  a given  time  as  In  Figure  4. 
The  amplitude  of  the  discontinuity  will  decrease  with  X. 


■V  Figure  (4) 

At  a given  point,  say  X >=  XQ  , the  variation  in  pressure  with 
time  is  indicated  In  Figure  5. 


completely  reasonable  assumption  if  the  period  is  very  much  greater 
than  the  mean  ilma  between  molecular  collisions.  On  the  other  hand, 
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the  change  from  b to  c Involves  an  Increase  In  entropy  as  calculated 
in  the  previous  section.  For  unit  cross  section  of  the  wave  the  in- 
stantaneous rate  at  which  entropy  increases  as  a wave  passes  Xn  can 
be  obtained  from  (9)  and  is  given  by 


and  hence, 


Xb 

XT 


— / 0 Co  ( ^ \ ) ■ 

_L  4J.  . 

- - ~ r d*  ‘ 


ct  £ 
dt 


WCL»  0 


til  s\ 

*" 


where  E is  the  mechanical  energy  of  vibration  per  wave  length  per 
unit  cross  section  and  p is  the  average  pressure.  In  writing  (10), 

SL 

the  velocity  of  the  wave  is  taken  an  c , where  c is  *he  sound  ve- 

o o 

locity  at  the  mean  temperature  of  the  gas.  This  can  be  justified  if 

it  is  noted  that  U is  just  equal  to  the  mean  value  of  u + c on  both 

sides  of  the  shock  to  within  first  order  terms.  * 

The  energy  E per  wave  length  per  unit  cross-  section  is 

2 

c.  — & ~ •/’Ot.wi.v  j 

& " — - p; • A ■ (12) 


With  the  assumption  about  the  shape  of  the  wave, 

'fa  fa-  & ( _ 

W — fly-  J xMx  — 


fait 


"“See,  e.  g.,  p.  159  of  Kef.  3. 
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Pace  12 


and 


/pa.  £ A 

1 3~  Y 


(13) 


Thus,  using  (11)  and  (13), 

I _ du% 

E d± 


— - -t.  A 


and  from  this 


i j.  ■£ 

e"  X* 


m j_ 

t'a  ’ 


In  terms  of  the  rate  of  change  of  & , 

j_  _ yy/  , 

b IT*  ~ y A 


(14) 


(IB) 


(10) 


If  (10)  Is  integrated  an  Interesting  relation  ensues,  namely, 
/ / V/-/  /"X, 


<2  -/ 


A 


(17) 


u. 


Discussion. 


of  wave  lengths, 


It  is  seen  from  (17)  that  a plot  of 

A -X® 


A 


-j-  against  the  number 
yields  a straight  line  whose  slope  is  wholly 


determined  by  / . It  is  to  be  noted,  as  one  would  expect  from  the 
nature  of  the  initial  assumptions,  that  the  attenuation  rate  is  independ- 
ent of  the  attenuatlve  process  or  the  constants  which  characterize  it. 
The  critical  assumption  is  that  the  wave  form  becomes  and  remains 


saw-tooth  in  character.  Thus,  any  attenuativ*  process  which  results 
in  the  presence  and  stability  of  the  saw-tootn  wave  form  will  lead  to 


process. 

It  Is  worth  emphasizing  that  the  only  entropy  change  which 
occurs  in  a wave  length  is  associated  with  the  shock  itself,  and  that 
no  entropy  change  occurs  on  the  gradually  sloping  part  of  each  wave. 
Consequently  for  waves  which  depart  slightly  from  the  assumed  saw- 
tooth form  (hut  which  nevertheless  possess  a shock  front  and  a grad- 
ually sloping  part)  of  complete  importance  lu  ^tdculating  the  rate  of 
energy  loss  is  the  magnitude  of  the  discontinuity.  The  nature  of  the 
gradually  sloping  part  of  the  wave  is  involved  only  In  the  determina- 
tion of  the  fraction  of  the  wave  energy  the  aforementioned  power  loss 
represents.  *. 


In  this  section  the  results,  equations  (1.5),  (16),  and  (17), 


will  be  discussed  in  the  light  of  earlier  related  work  by  Dumond, 
Cohen,  Panofsky  and  Deeds*,  and  Fa^ . 

4 

Dumond  et  al.  derive  the  attenuation  rate  for  the  N wave 


CO 
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where 


They  find  that 

iii- 
e d. 7 


^3  __ 


-s.  I + L * 


ItL 

^y 


A 

A 


(18) 


This  result  could  also  be  obtained  by  the  same  methods 
used  here.  Applying  (11)  to  the  propagation  of  the  two  shock  fronts 
of  the  N wave,  one  obtains  by  adding  the  two  equal  loss  rates, 


if  - 


and  since  E = A^-sy  > 


a* 

i/1 

it  is  seen  that  Eq.  (18)  follows. 


As  pointed  out  by  Dumond  et  al. , the  propagation  of  the  N 

wave  is  accompanied  by  a lengthening  of  the1  wave.  That  this  is  so 

is  cloar  from  the  fact  that  tire  velocity,  U,  of  the  forward  shock  is 

given  by  the  average  of  u + c on  the  two  sides  of  the  shock.  It  is 

y f i 

thus  equal  to  c9  ■+■  ^ • Similarly,  the  velocity  of  the  rear 

shock  is  — y 'LL  u.  , Thus, 

d A / //  (j  — y+i  a 

Tt  ” — u ~ z-i  ‘ ' 

where  the  relation  u = ^eC.  has  been  used.  For  the  re- 

pealed shock  wave-,  on  the  other  hand,  the  velocity  of  consecutive  fronts 
is  the  same  to  within  terms  of  the  first  order  in  & , and  consequently 
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the  wave  length  is  constant. 

It  is  worth  while  showing  that  the  alternate  hydrodynamic 
method  Dumond  et  al.  use  to  derive  Eq.  (18)  may  be  very  nicely  ap- 
plied to  the  repeated  shock  and  yields  Equations  (15)  and  (10).  Con- 
sider a linear  pressure  gradient  in  a sound  wave  represented  by  AB. 
Then,  if  0 represents  the  point  at  which  the  particle  velocity  Is  zero, 


ttludt 


A P"  i . 7 . 

every  given  point  on  the  line  AB  propagates  away  from  this  point  with 
a velocity  j where  _u  is  the  particle  velocity  at  the  given 

point  in  the  wave.  (The  point  0 propagates  with  the  velocity  of  sound, 
cQ  . ) This  is  a result  of  solving  the  fundamental  hydrodynamicai 
equations.  Thus,  after  a time  interval,  dt,  the  pressure  curve  would 
be  described  with  respect  to  the  displaced  coordinate,  + C0<iX  . 
However,  if  there  are  shock  fronts  DE  and  FO,  so  that  ODE  represents 


one  wave  of  a repeated  shock  of  the  typ>  being  discussed,  then  DE  will 
propagate  vrtth  a velocity  which  is  the  mean  value  of  c + u on  both 
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sides  of  the  shock  which  la  equal  to  cq  . Thus,  a wave  which  satis- 
fies both  requirements  Is  one  whose  typical  wave  length  is  represented  by 
LKIIJ.  (Ho  is  the  half  height  of  the  discontinuity. ) Nov;  DJ  - /f°a.^  ) 


and  DU  - - cL  (fix  . From  the  geometry. 


d-  ('Kx) 


11L  licit 

Whence,  using  ^ ^ 

Fq  C*  O 

(10)  is  derived. 


**  m ljt±  n,,d  Celt  =■  d* 

° Po 


Z * 

Fay's  result,  Eq.  (14)  , can  be  used  to  derive  (16)  of  this 
paper.  To  show  this,  his  equation  is  recast  in  terms  of  the  variables 
used  here. 

ft -.bo  jlL.J..2 iL)  y ^ ^ ^ ) 

~ It!  «3  /?  Co 


— - * * it#  ' 

= kinematic  coefficient  of  viscosity  } 

* 0 


ui  = 7ir 


X > 


<U  o - constant. 


As  an  example,  determine  the  value  of  °<c  , C.X  = Xc  ) , If  the 


Note  that  the  factor  8 in  Eq,  (14)  of  Ref.  2 should  be  2. 
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fundamental  of  the  wave  la  to  have  a maximum  value  of  ^tlil  - ao  5 , 


UJ  z £ oo 


w 

77  , 7r  / H J S^//>a  - d./3  , Cc  =:  3,  5 ‘ V 


It  la  found  that  ^ ^ wwi  ^ * V * / ^ 


- £ 


^ 5“ 

Thus,  for  Fv\  < | £ , the  hyperbolic  sine  can  be  replaced  by  its 

argument  and  the  amplitude  of  the  harmonics  is  inversely  proportional 
to  the  order  of  the  harmonic,  all  harmonics  being  in  phase.  More- 
over,  since  U.  = O.^U'  10  j « X < oifl  for  X < 4 meters.  With 
conventional  measuring  instruments  the  wave  will  appear  to  retain  its 
saw-tooth  character  for  many  wave  lengths  of  travel.  For  example, 
if  a microphone's  response  drops  off  seriously  above  10  kc  (the  100th 
harmonic  of  the  100  c.p.  s.  wave),  the  calculation  indicates  that  the 
distance  of  propagation  before  a significant  departure  from  saw-tooth 
character  would  be  observed  is  approximately  15,000  meters  (assuming, 
of  course,  plane  wave  propagation). 

To  continue  with  the  original  purpose,  suppose 

\ o.o3 

V f)  4 A >'M( 

Then  till  ^ A-  K tH  1 ( — " JEI  , 

4*  c IT 

-w  t 

where  o is  giver,  by  Fq,  (8). 
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i w*i  A'i.,(  uALir.ia.TA 

TT  _ ttL  *)+•** 

"T  ” 2 y a 

3 /*  3- 

' o ^ c 

Hi)  = , 

which  checks  Eq.  (18). 

It  is  perhaps  worth  emphasizing  a point  which  has  been  made 
by  Fay  — if  one  adopts  the  approach,  which  most  people  in  acoustics 
think  of  naturally,  of  reducing  a wave  to  its  Fourier  components,  then 
the  attenuation  of  a saw-tooth  wave  occurs  principally  in  its  high  fre- 
quency components.  The  procedure  by  which  a saw-tooth  wave  main- 
tains its  form  is  by  the  low  frequency  components  feeding  energy  into 
the  higher  harmonics  and  the  higher  harmonics  losing  vibrational  en- 
ergy through  the  action  of  viscosity.  Thus,  in  Eq.  (10)  the  term 
oac*)*,  in  the  numerator  is  associated  with  attenuation,  and  a large 
value  of  (e.  g. , o'  » If  in  the  example  taken)  is  evidence 

of  a large  rate  of  attenuation. 

Another  7/ay  of  looking  at  this  is  to  calculate  the  rate  of 
energy  loss  of  a saw-tooth  wave,  which  maintains  its  shape,  by  as- 
suming that  the  total  rote  of  energy  loss  is  the  sum  of  those  due  to 
the  harmonics,  assuming  each  harmonic  is  attenuated  at  a fractional 
rate  which  is  given  by  Infinitesimal  amplitude  acoustics.  Let  F^  be 
the  vibrational  energy  density  associated  with  the  nth  harmonic,  and 


be  the  total  energy  density.  Then 
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With  our  assumption, 
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where  <X,is  the  value  of  <X  (defined  below  Eq.  (19)  ) for  the  funda- 


mental. 
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The  summation  in  the  numerator  does  not  converge.  On  the  other 
hand,  if  the  assumption  is  made  that  the  harmonics  of  order  greater 


than  on  --  N are  absent,  then, 
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monies  N which  is  necessary  to  give  the  rate  of  attenuation  given  by 


Eq.  (To). 
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For  a fundamental  of  100  c.p.  s.  and  j = 0.03, 
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Since  there  are  more  than  ten  times  as  many  harmonics  between  1 and 
13  megacycles  as  below  1 megacycle,  til's  calculation  would  indicate 
that  most  of  the  energy  lost  by  viscous  action  occurs  for  components 
above  1 megacycle.  While  the  details  of  this  calculation  cannot  be 
correct,  this  last  conclusion  is  correct  with  better  than  order  of  mag- 
nitude accuracy. 
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